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Abstract 

It is analyzed the hydrodynamics of carriers (charge and heat motion) and phonons (heat mo¬ 
tion) in semiconductors in the presence of constant electric fields. This is done in terms of a so- 
called Higher-Order Generalized Hydrodynamics (HOGH), also referred to as Mesoscopic Hydro- 
Thermodynamics (MHT), that is, covering phenomena involving motions displaying variations 
short in space and fast in time and being arbitrarily removed from eqnilibrium, as it is the case 
in modern electronic devices. The particular case of a MHT of order 1 is described, covering wire 
samples from macro to nano sizes. Electric and thermal conductivities are obtained. As the size 
decreases towards the nanometric scale, the MHT of order 1 produces results that in some cases 
greatly differ from those of the usual hydro-thermodynamics. The so-called Maxwell times associ¬ 
ated to the different fluxes present in MHT are evidenced and analyzed; they have a quite relevant 
role in determining the characteristics of the motion. 
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I. INTRODUCTION 


The modern advanced technologies, and its resulting end use for improved and novel 
products, create a stress on the basic sciences of Physics and Chemistry. This is a result 
of trying to mantain a balance in the triade ST&I (Science, Technology and Innovation) 
[1], Particular questions involve, for example, the dissipation of energy and heat transport 
in devices under high-levels of excitation, namely, working in far-removed-from equilibrium 
conditions and eventually involving ultrafast relaxation and transport processes, as well as 
spatial motion in nanometric scales. Another important aspect is the one of fluids under 
flow present in certain production processes (e.g., in food engineering, petrochemistry, etc...) 
whose performance depends on their hydrodynamic properties [2]. Moreover it can be men¬ 
tioned the question of hgure of merit in thermoelectric devices, that is, relating currents of 
charges and of heat, particularly in the nanometric scale [3]. 

It has been noticed that one of the complicated problems of the nonequilibrium theory 
of transport processes in dense gases and liquids is the fact that their kinetics and hydro¬ 
dynamics are intimately coupled, and must be treated simultaneously (e.g., see Refs. [4-6]). 
Along the last decades Hydrodynamics has been extensively treated resorting to the so-called 
Nonequilibrium Molecular Dynamics (NMD for short). NMD is a computational method 
created for modelling physical systems at the microscopic level, being a good technique to 
study the molecular behavior of several physical processes [7,8]. On the other hand, another 
very satisfactory approach to deal with hydrodynamics within an ample scope of nonequi- 
librium conditions consists in the kinetic theory based on the Non-Equilibrium Statistical 
Ensemble Eormalism (NESEF for short) [9-14]. NESEF is a powerful formalism that provides 
an elegant, practical, and physically clear picture for describing irreversible processes, as for 
example in semiconductors far-from equilibrium [15-17]. NESEF provides a way to go be¬ 
yond standard (or classical) Onsagerian hydrodynamics which involves restrictions, namely, 
local equilibrium; linear relations between fluxes and thermodynamic forces (meaning weak 
amplitudes in the motion) with Onsager’s symmetry laws holding; near homogeneous and 
static movement (i.e., involving only smooth variation in space and time); and weak and 
rapidly regressing fluctuations [18,19]. Hence, more advanced approaches are required to 
lift these restrictions. In phenomenological theories this corresponds to go from classical 
irreversible thermodynamics to extended irreversible thermodynamics [20-22]. This is what 
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has been called generalized hydrodynamics, a question extensively debated for decades by 
the Statistical Mechanics community. Several approaches have been used, and a descrip¬ 
tion can be consulted in Chapter 6 of the classical book on the subject by Boon and Yip 
[23]. Introduction of nonlocal effects for describing motions with influence of ever decreas¬ 
ing wavelengths, going towards the very short limit, has been done in terms of expansions 
in increasing powers of the wavenumber, which consists in what is sometimes referred to 
as “Higher-Order Generalized Hydrodynamics” (HOGH for short), also dubbed Mesoscopic 
Hydro-Thermodynamics (MHT for short) [24]. 

Within the scope of Mesoscopic Hydro-Thermodynamics we consider here the question of 
transport of charge and of heat in n-doped polar semiconductors in the presence of electric 
helds. The hierarchy of equations of evolution for the density and energy density of carriers 
and of energy density of phonons and together with those for their fluxes of all orders, 
are obtained in the framework of the nonlinear quantum kinetic theory that is based on 
NESEF [11-13,25,26]. The electrical and thermal conductivities in such nonequilibrium 
thermodynamic state and within a MHT of order 1, that is, a description reduced to include 
the densities and their first fluxes, are derived, and the influence of the order of the HOGH 
(contracted description in terms of the densities and a reduced number of higher-order 
fluxes) and of the sample size (macro to nano) are discussed. The so-called Maxwell times 
[27,28] are characterized and analyzed, with some numerical calculations and hgures being 
presented. Maxwell times are of fundamental relevance for establishing the order of the 
contracted description of MHT to be used, and of large influence on the behavior of transport 
properties at short nanoscales. 


II. THEORETICAL BACKGROUND 

The construction of a Mesoscopic Hydro-Thermodynamics for the description of the move¬ 
ment of matter and energy in fluids under nonequilibrium thermodynamic conditions and 
at the classical mechanical level based on a generalized moments approach method to the 
solution of a NESEF-based generalized Boltzmann equation [29], is described elsewhere 
[30,31]. 

We consider here MHT at the quantum mechanical level, for dealing with a system of 
carriers and phonons in n-doped polar semiconductors in the presence of electric fields (up to 
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100 kV/cm) which drive the system away from equilibrium. Moreover, the system is taken 
to be in contact with an external thermostat at temperature Tq. 

The system is characterized at the microscopic level by the Hamiltonian 

H = He + Hp + Hee + Hep + Han + He£ + HpR , (1) 

consisting of the Hamiltonians of the free electrons and free phonons, respectively 

He = Yl CkC^Ck , ( 2 ) 

k 

Hp = Y^ + 1/2), (3) 

q7 

where Ck is the electrons’ conduction-band energy (spin index has been ignored), cUq.^ is the 
phonon frequency dispersion relation with 7 indicating the branch LO, to, la, ta, and k 
and q are wave vectors running over the Brillouin zone. The electron-electron interaction is 

Hee= H(kik 2 k 3 k 4 ) 4 ^ 42 Ck 3 Ck 4 , (4) 

kik2k3k4 

(ki+k2=k3+k4) 

and for the electron-phonon interaction we have 

~ I'kq7^q7*^k+q'^k + H-C. . (5) 

kq 

where as noticed 7 indicates the phonon branch, and a the type of interaction (deformation 
potential, Frohlich-polar with LO phonons, piezoelectric with LA phonons) with coupling 
strength C. The electron-electric field interaction is given by 



with an electric held E of intensity S applied in, say, ; 2 -direction and we write for the 
anharmonic interaction 

Han = 'y ^ Tfkq7&q7&k4.q^/&Lk7'' + H-C- , (7) 

kq 

ll'l" 

(where we have neglected nonlinear contributions), with Mk^q^ accounting for the coupling 
strength. 
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For the description of the macroscopic nonequilibrium thermodynamic state we resort, as 
already noticed, to the use of NESEF. The statistical approach NESEF requires hrst of all 
to specify the basic dynamical variables used to characterized the non-equilibrium ensemble 
[9,11-13,25,32], A priori, when the system is driven away from equilibrium, it is necessary 
to include all observables of the system, which leads to the introduction of many-particle 
dynamical operators [33,34], in the present case of single electrons in Bloch conduction band 
and single phonons, it sufhces to introduce only the single particle dynamical operator, 
namely 

hk ) hkQ *^k-|-Q/2*^k—Q/2 ) (8) 

with Q 7 ^ 0, and spin index ignored, for the carriers, and 

f'qi = ^q7^q7 i ^qQ7 “ ^q+1 Q,y^q-iQ7 ’ (9) 

with Q 7 ^ 0 for the phonons. 

Dynamical operators of order two and higher (in the BBGKY hierarchy [34]) do not 
contribute because correlations and higher-order variances are absent in the mean-held ap¬ 
proximation for the carriers and the harmonic approximation for the lattice vibrations. 
Moreover, since phonons are bosons, it would be necessary also to include the annihilation 
and creation operators 6 q and because their eingenstates are the coherent states [35], 
and also the pair operators &k&k', ^iL^k' because the number of quasi-particles is not hxed 
[36]. However, we disregard them because are of no relevance for the problem at hands. In 
Appendix A we describe the corresponding non-equilibrium statistical operator (Cf. Eqs. 
(Al) and (A 2 )). 

Operators hk and correspond to the occupation number operator describing a homoge¬ 
neous population and those with Q 7 ^ 0, account for changes in space of the non-equilibrium 
distribution functions. 

The average, over the non-equilibrium ensemble, of the microdynamical variables in the 
sets of Eqs. ( 8 ) and (9) provide the variables which characterize the non-equilibrium macro¬ 
scopic state of the system. Let us call them 

{'^k(^))'^kQ(^)) ^q7(^)) ^qQ7(^)} i (10) 

where Q 7 ^ 0 and n\^{t) = Tr {r^k^?£(t)}, etc..., that is, the average over the non-equilibrium 
ensemble according to the formalism in Appendix A where we have introduced the non- 
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equilibrium thermodynamic state variables said conjugated to those above, namely [cf. Eq. 
(A.3)] 

. ( 11 ) 

Going over to direct space we introduce the space and crystal-momentum dependent 
distribution functions 


^5^nkQ(t)e*^-" = /k(r,t), 

(12) 

= ^q(rC), 

Q 

(13) 


where Vcei is the volume of the unit cell, and the phonon branch under index 7 is being 
implicit from here on. 

In terms of this microscopic (quantum mechanical) and macroscopic (nonequilibrium 
thermodynamic) description of the system, we proceed to present the evolution equations 
of the basic variables in the set of Eq. (10). 

III. EVOLUTION OF THE NONEQUILIBRIUM THERMODYNAMIC STATE 

Calling, in a compact and generic form, Pj and Qj{t) the dynamical variables and the 
corresponding thermodynamic variables in the sets of Eqs. (8) and (9) and (10) respectively, 
the evolution equations for the variables Qj {t) describing the evolution of the nonequilibrium 
thermodynamic state of the system are 

^<5j(^) = X 0b} = Tr H]ge{t) x , (14) 

that is, the average over the nonequilibrium ensemble, characterized by the statistical op¬ 
erator Qsit) of Appendix A (Cf. Eqs. (Al) to (A3)), of the Heisenberg equation for the 
corresponding dynamical variable Pj; qb is the distribution of the surrounding medium 
assumed in equilibrium at temperature Tq. 

Direct calculation of the right-hand-side in Eq. (14) is extremely difficult and then 
it is necessary to resort to the introduction of a more practical non-linear kinetic theory 
[11,13,25,26] briefly described in Appendix B, which is applied using an approximation 
consisting in retaining only the collision integrals of second order in the interaction strengths 
(Markovian approximation [26,29,32]). In reciprocal space it follows that (see Appendix B) 

= -^(^Q ■ Vkek)nkQ(t) -h , (15) 
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(16) 


= W ■ + <q(*) + <S(«) + (*). 

where the seven collision integrals J’s are given in generic form in Eqs. (B.3) to (B.6) in 
Appendix B. In direct space, after using Eqs. (12) and (13) and for the different J’s that 


Jk(r,t) 


1 


Q 


; A(r,t) 


1 

hceZ 


Q 


we do have that 

J/klr, t) + 1 Vkfk ■ V,A(r, t) = Jf "”“(r, (), 

where: 

Jf *““(r,«) = Jf"(r,«) + Jr(r,«) + Jf*(r, t ), 


Jj"””‘(r.i) = Jr(rJ) + Jr(rJ) + JVM . 


qj 


qj 


qj 


( 17 ) 

(18) 

(19) 

( 20 ) 
( 21 ) 


where J^'^(r,t) accounts for the effect of the electron-phonon interaction, J^'®(r, t) of the 
internal interaction (electron-electron), J^'^(r,t) the electron-sources interactions, J^^{r,t) 
of the phonon-electron interaction, Jq)) (r,t) of the anharmonic interaction and of 

the phonon-sources interactions. 

We consider now the Mesoscopic Hydro-Thermodynamic of the system, which consists 
into introducing the densities of (quasi)particles and of the energy and their fluxes of all 
order, namely, for the electrons. 


{n,(r,«),I„,(r,i),.../W(r,«)...}, 


( 22 ) 


which we call MHT-carriers’ family n, and 

{he(r,t),I/,^(r,t),.../J2(r,f)...}, 
the MHT-carriers’ family h, where 

^e(r,t) = ^/k(r,t) , 

k 

Ine(r,t) = ^ Vkek/k(r,t), 

k 

k 


(23) 


(24) 

(25) 

(26) 
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with £ = 2, 3,... and 


“fk = ^[^kek: .. .^ times...; VkCk] , (27) 

is a rank-t' tensor involving £-times the tensorial internal product of the group velocity 
(l/h)Vkek (in a effective mass approximation Ck = h^k‘^/2m*, and then VkCk = Kk/ml] m* 
is the effective mass of the electrons at the center of the conduction Bloch band in polar 
semiconductors), and 

^e(r,t) = ^ek/k(r,t) , (28) 

k 

I/ie(r,t) = ^ekVkek/k(r,t), (29) 

k 

4! (!■> ^) = ek44/k(r, t). (30) 

k 

On the other hand, for the phonons we do have 

{np(r, t), I„p(r, t),..., /Jg(r, f),...} , (31) 

the MHT-phonons’ family n, and 

{hp{r,t), I,,^(r, t),..., 4^j(r,t),...} , (32) 

the MHT-phonons’ family h, where 


np{^,t) = ^z/q(r,t), 
q 

(33) 

Ip(r, t) = ^ fhUqVqWq Z/q(r, t) , 
q 

(34) 

4^’ (r, t) = Yl ^q41q ’ 

(35) 


q 

for the MHT-phonons’ n-family, and 

hp(r,t) = ^ ha;qZ/q(r,t), (36) 

q 

Iftp (r, t) = Y^ hWqVqWq Z/q(r, t) , (37) 

q 
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(38) 



for the MHT-phonons’ h-family, and where 

“S,q = [VqWq I ... £ timeS . . . ; VqWq] , (39) 

and we recall that the phonon branch index 7 is implicit; VqWq is the group velocity of 
phonons in mode q. 

The evolution equations which describe the hydrodynamic motion in MHT are: 


^4!!(r,t) = 5^MW(k)^/k(r,t), 

k 

(40) 

= 5^ekM|fl(k)^/k(r,t), 

k 

(41) 

= 5^Mjl(q)^^^q(r,t) , 

q 

(42) 

= 5^^qWpl(q)^'^q(r,t), 

(43) 


q 


where, we recall, £ = 0 stands for the densities, £ = 1 for their first (vectorial) fluxes, 
£ = 2,3,... for the higher-order tensorial fluxes. 

This set of equations is practically intractable, requiring to look in each case on how to 
hnd the best description using the smallest possible number of variables. In other words 
to introduce an appropriate - for each case - contraction of description: this contraction 
implies in retaining the information considered as relevant for the problem in hands, and to 
disregard irrelevant information [37]. 

Elsewhere [38] it has been discussed the question of the contraction of description (re¬ 
duction of the dimensions of the nonequilibrium thermodynamic space of states), where a 
criterion for justifying the different levels of contraction is derived: It depends on the range of 
wavelengths and frequencies which are relevant for the characterization, in terms of normal 
modes, of the hydro-thermo dynamic motion in the nonequilibrium open system. It can be 
shown that the truncation criterion rests on the eharacteristies of the hydrodynamie motion 
that develops under the given experimental procedure. 

Since inclusion of higher and higher-order fluxes implies in describing a motion involv¬ 
ing increasing Knudsen numbers per hydrodynamic mode (that is, governed by smaller and 
smaller wavelengths - larger and larger wavenumbers - accompanied by higher and higher 
frequencies). In a qualitative manner, we can say that, as a general “thumb rule”, the crite¬ 
rion indicates that a more and more restricted eontraction can be used when larger and larger 
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are the prevalent wavelengths in the motion. Therefore, in simpler words, when the motion 
becomes more and more smooth in space and time, the more reduced can be the dimen¬ 
sion of the basic macrovariables space to be used for the description of the nonequilibrium 
thermodynamic state of the system. It can be conjectured a general contraction criterion, 
namely, a contraction of order r (meaning keeping the densities and their fluxes up to order 
r), once we can show that in the spectrum of wavelengths, which characterizes the motion, 
predominate those larger than a “frontier” one, >^lr,r+l) = V^0r9r+1 where v is of the order of 
the thermal velocity and 9r and 9r+i the corresponding Maxwell times, see next, associated 
to the r and r -|-1 order fluxes. We shall try next to illustrate the matter using a contraction 
of order 1, that is, a first-order extension of standard Onsagerian hydrodynamics. 


IV. MHT OF ORDER 1 OF CARRIERS AND PHONONS 

We consider the contracted MHT of order 1 (that is keeping only the densities and their 
first fluxes, implying in smooth movement in space and slow in time, but beyond the range 
in standard hydrodynamics) in the already described system of carriers and phonons in a 
n-doped polar semiconductor in the presence of an electric held. Hence, the basic sets of 
dynamical variables are 



|he(r),i„^(r),he(r),ife^(r)| 

(44) 

for the carriers, and 

{hp(r),i„^(r),hp(r),i,,^(r)| , 

(45) 

for the phonons, or in 

reciprocal space 



|he(Q),i„,(Q),he(Q),i/ie(Q)} , 

(46) 

and 

|hp(Q),i„^(Q),hp(Q),i/ip(Q)| . 

(47) 
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The associated auxiliary statistical operator g (see Appendix A) is then 
g{t, 0 ) = exp I - 0(t) 

~ ^)^e(Q) + F„^(Q, t) ■ ine(Q; t) + 

Q 

+TX^(Q, t)he{Q„ t) + F/i^(Q, t) ■ ihe(Q) t) + 

+-Fnp(Q,t)hp(Q) + F„p(Q, t) ■ i„p(Q) + 

+-Fftp(Q,t)hp(Q) + Fftp(Q,t) ■ i/jp(Q)]| , (48) 

introducing the set of nonequilibrium thermodynamic variables 


where Q = 0 refers to the homogeneous (or global) state-thermodynamic variables, and 
Q 7 ^ 0 to the inhomogeneous contributions. As usually done, we write 




F,Mt) ^ ^ -[kBT:{t)]-\ 


(50) 


introducing the carriers quasitemperature T*(t), 


Fn,(0,t) = /3e{t)\e{t) , 


(51) 


with Ve(t) being the drift velocity, and 


Fn,(0,t) = /3e{t)fil{t) , 


(52) 


introducing the quasi-chemical potential /i*. 

On the other hand, for the phonons we do have that Fnp{0,t) = 0 (number of phonons 
not conserved), F„p(0,t) (no close current circuit present), and we write 

Fh^i0,t)hp{0) = , (53) 

q7 

with 

/?„(() = |fer*,(()]-‘, (54) 

introducing the phonon quasitemperature per mode in each branchT^^ (t). 

The set of evolution equations, see Appendix B, for the electrons results in that 

^ne(Q, t)=iQ- In, (Q, t) + <he(Q, t) , (55) 
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d 

dt 


= tQ- lS{Q,t) 

^ ^ (^k+q*^k+q VkCk)^kq^o(^) (^e,kQ(^) ^e,k,— q(^)) 

kq7Q: 

+eEne(Q,t), 


( 56 ) 


dt 


he{Q,t) = iQ ■ lh,iQ,t) 


^ ^ (^k+q Ck)^kq^a(^) ("^ejkQ (^) “1“ '^e,k,— q(^)) 

kq 7 Q: 

E 

+e— ■ Ie(Q 7 )^e(Q 7 ) , 

ml 


(57) 




*Q'iS(Q.«) 

^ ^ (^k^k^k ^k+qVk+q^k+q)^kq^Q((^)'^kQ(^) 

kq7Q; 

+lPl(q,«).eE, 


(58) 


where 


O-TT 

^kq7,aW = 7-|Ckq^P{[(^^q7W + l)/k+qW + ^^q7W(l-/k+qW)] 

X<^(Ck+q Ck hcjq-y) 

+ [(z/q 7 (t) + 1)(1 — /k+q(t)) + r'q 7 (t)/k+q(t)] 

x5(ek+q - ek + hcj q7)} > 


(59) 


with the presence of the homogeneous distributions 

m = ^ 


e^e(t)[ek-/7hb-Ve(t)'k] _j_ ]_ ’ 


( 60 ) 


resembling a kind of shifted instantaneous Fermi-Dirac distribution, which in the nondegen¬ 
erate limit becomes 


fUt) = 4nW ( ] X 

2m* 


(61) 


i.e., a shifted Maxwell-Boltzmann-like distribution where n is the density of carriers, and 


t+\ ^ 


(62) 
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which has the form of a Bose-Einstein like distribution at zero quasi-chemical potential. 
Equations (59) and (61) are a result of the calculation of 

/k(t) = Tr{ 4 ckp(t, 0 )} , and z/q^(t) = Tr{&]j^ 6 q^p(f, 0 )} , 

and the use of Eqs. (50), (51), (52) and (54). 

The scattering integral <he in Eq. (55) accounts for local effects due to the presence 
of impurities, imperfections (dislocations, stacking faults, etc.), the imperfections in the 
end contacts, and geometry and boundary influences; the lateral walls are rugous (of a 
fractal-on-average character [39]) leading to inhomogeneous scattering of the carriers. The 
integration of ‘he(r, t) over the volume of the sample is null since the total number of carriers 
is constant. We recall that Ckq 7 matrix element of the electron-phonon interaction 

[cf. Eq. (5)]; it may be noticed that in polar semiconductors Frohlich-polar interaction (7 = 
LO, a = Frohlich interaction) is by far the relevant one producing rates of change orders 
of magnitude greater than those associated to the other interactions [40,41]. Moreover, we 
have neglected the contribution of the plasma states via Coulomb interaction. 

For the phonons we do have, 

J^np{Q,t) = iQ ■ In,(Q,t) - ^ 5^[rq^Q/2,7a(^) + 

-5 + rr-Q/ 2 ,,„(*)]!'«,(*). (63) 

q'ya 

jfKW) = *QvP'(Q.*) 

q'ya 

-5 E|rr+Q/2,7<.(‘) + r“”Q/2,7„(«)lVqU^q,-2qQ7(«) , (64) 

q'ya 


^hp(Q,t) = iQ-I;,p(Q,f) 

5Z[rq+Q/2,7«(^) + rq!Q/2,7„(^)]^q7^qQ7(^) 

qya 



q 


(65) 
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~2 ^q7^q^q7[^q+Q/2,7o(^) ^q-Q/2,7«(^)]^qQ7(^) 

q7Q; 

-^5^^q7'^q '^q7[^q+Q/2,7«(^) rq_Q/2,7a(^)]t'qQ7(^) , (66) 

q7Q; 

where 

KU*) = EiQf-,P{/k'+,(«)ii-/k'Wi-/k-(«)ii-/k'+,(«)]} X 

k' 

^(ek'- Ck -/hk;q7), (67) 


with dimension of inverse of time, accounts for the rate of transfer (energy and momentum) 
from the hot carriers, and 

TT 


q7Q; 


it) = 




ElK7q'7'l^(l + Z^q/ 7 / l^d+qf'-y )Si ^q+q'7 + l^q' 7 ' I^q 7 )+rq, ( 68 ) 


with the first contribution on the right being the explicit expression for the inverse of the 
relaxation time due to anharmonic interactions, and F stands for, in a Mathiessen-like 
rule, the sum of the reciprocals of the relaxation times associated to the interaction with 
impurities, imperfections, stacking faults, as well as effects of (rugous) boundary conditions 
and contacts with other subsystems and sources. 

Evidently, the set of Eqs. (63) to (66) is not a closed one, once the right-hand sides are 
not given in terms of the proper basic hydrodynamic variables. Hence, we must proceed to 
introduce a closure condition, what is done resorting to Heims-Jaynes perturbation procedure 
for averages [42]. This is described in Appendix C, and in a first-order linear approach in 
Heims-Jaynes procedure, the fundamental set of hydrodynamic equations in MHT of order 
1 are for the carriers 

^ne(Q, t)=iQ- In, (Q, t) + d>e(Q, t ), (69) 


4 i„.(Q 7 ) 


HgiQne(Q, t) + (Q, t) - Oj^In^Q, t) 

eE 

+&34eI/ie(Q) f) J-, 

m ^ 


(70) 




iQ ■ Ih,(Q,t) + b2iene(Q,t) - 7/he(Q,t) 
-heE ■In,(Q,t)ne, 


( 71 ) 
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( 72 ) 


where 


with 


d 

dt 


^he{Q,t) 


C|giQne(Q,t) + C2e^Qhe(Q,t) + 643eIne(Q,^) 
-ej^lhAQ, t) + Bf}ne{Q, t)E + t)E, 


Bfht) = 5^[Vkek; Vkek]&ie(k,t), (73) 

k 

B?}{t) = 5^ek[Vkek; Vkek] 62 e(k,t), (74) 

k 

cSit) = 5^ek[Vkek; Vkek]6ie(k,t), (75) 

k 

C?ht) = 5^[Vkek: Vkek] 62 e(k,t), (76) 

k 

6ie(k,t) = [A 12e it)r[A 22e{t) — ek^l2e(t)]/k(t)[l — /k(t)] , (77) 

&2e(k,t) = [A 12e lle(^)ek — ^21e (t)]/k(t)[l-/k(t)], (78) 

7liie(t) = 5^/k(t)[l-/k(t)], (79) 

k 

^12e(t) = A21e(t) = ^ek/k(t)[l — /k(t)] , (80) 

k 

7l22e(t) = 5^(ek)Vk(t)[l-/k(t)], (81) 

q 

^12eit) = 24iie(t)A22e(t) — Ai2e(t)^216(t) . (82) 


with /k(t) of Eq. (60), and we recall that [... : ...] stands for inner tensorial product of two 
vectors producing a rank -2 tensor. 

It can be noticed that these expressions can be greatly simplihed if we disregard the 
contribution of the drift velocity Vg in the distribution /k(t) (the kinetic drift energy is 
smaller than the thermal energy for any intensity of the electric held [43,44], and then, 
because of the spherical symmetry in the expressions for the tensorial kinetic coefficients, 
they become scalars. 

Moreover, in Eqs. (70), (71) and (72) are present generalizations of the so-called Maxwell 
time [27,28], 9h,, ^ 4 ^, given by 

[K. (^)]"^ = 5 ^ ■ k7lkq/k(l - /k) , (83) 

f I i 
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( 84 ) 


where 


[^/ie(^)] ~ ^ ^ ^kq^2e(Ck+q ^k) ; 

k,q 

(^)] ~ ^ ^ ^kqb4e(ek^k€k ^k+qVk+qCk+q) ) 

k,q 


^^kq ^ l^kql {[(^q “1“ l)/kH-q “1“ ^q(l /k+q)]<^(Ck+q ^k fi^q) 

[(^q “1“ 1)(1 /k+q) “1“ ^q/k+q]<^(ck+q ^k “1“ hjJq) j 


(85) 


( 86 ) 


&2e is given in Eq. (78) and 

b4e = [A34] ^[AsaekVkCk — A34Vkek]/k(l — fk), (87) 

A34 = A33A44 — A34A43 , (88) 

^33 = ^ |VkekP/k(l -/k), (89) 

k 

A44 = 5^(ek)'|Vkekr/k(l - fk), (90) 

k 

A34 = ^ ek|VkekP/k(l —/k) = A43 . (91) 

k 


Neglecting the dependence on time of all the different coefficients (i.e., taken them as 
weakly dependent on time), going over direct space, the basic equations of the MHT of 
order 1 of the carriers in doped semiconductors are 

^ne(r,t) + V ■I„^(r,t) = <he(r,t), (92) 


m 


Ine(r,t) 


B[^J(t)Vne(r,t) - 0jJ^In,(r,t) + b34(t)Ih, (r,t) + 
+41Wl/4e(r,t) + ^ne(r,t) , 


(93) 


m 


he{r,t) + V ■ Ih,(r,t) 


-9^^he{r, t) + h2iene{r, t) + 
+eE ■ I„^(r,t)ne(r,t), 


(94) 




I/4e(r,^) 


-CgVA.(r, t) - 97n»,(r, t) - ClfVn.(r, t) + 

+443e(<)In.(r,«) + i) + -BfjE/i,(r,«) . 


(95) 
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In Eq. (92), on the left-hand side is present the barycentric time differentiation (the 
conservation part), and on the right the source of local variations due to the presence of 
impurities, boundaries, etc. (we recall that the integration in space of it is null because the 
conservation in the number of charges). Equation (93) is on the right composed of a hrst 
term of a diffusive character, followed by Maxwell contribution, the third contribution is a 
cross term associated to thermo-striction effects, and the last one accounts for the effect of 
the presence of the electric held creating the electric current. 

In Eq. (94), the left hand side represents the conserving part of the energy, and on the 
right we hrst hud Maxwell contribution which is followed by a cross term associated to 
thermo-electric ehects, and a contribution due to the presence of the electric held. The last 
one is the local production of Joule heat. 

In Eq. (95) various terms contribute on the right: the hrst is of a dihusive character, 
followed by Maxwell contribution. The third and fourth terms are cross terms associated to 
thermo-striction ehects, and the last two are contribution due to the presence of the electric 
held. 

We consider next the associated hydrodynamic modes. 

V. THE HYDRODYNAMIC MODES IN MHT [1] 

For the purpose of obtaining the hydrodynamic modes of the carriers in the MHT of 
order 1, we consider Eqs. (69) to (72), but introducing the simplihcations of neglecting the 
source <he in Eq. (69), i.e., disregarding the ehect of impurities and imperfections. In the 
last term in Eq. (71) we take for ne(r, t) only the relevant constant uniform contribution nP, 
i.e., the doping concentration, and we take the second-rank tensors, Ci^ and , 

as scalars, all of this to have manageable equations for just to better visualize the physical 
characteristics of the hydrodynamic motion. 

Transforming Fourier in time Eqs. (69) to (72), we are left with a set of linear algebraic 
equations (in Q-u space) whose secular determinant is 


ioj 

—iQ 

0 

0 


iu + 

in 

-i?2*Q 

CO 

1 

—b2i 

— en°E 

iu + e-^ 

— fQ 

-CiiQ - BiE 

~&43 

-C2fQ - B2E 

lu + djl 
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The complete set of hydrodynamic modes are the solutions, say a;i, 2 , 3 , 4 , of a fourth-order 
algebraic equation, which follows after making this determinant equal to zero, which we 
omit to write down explicitly. We consider now a situation when thermo-electric effects can 
be neglected, and then the movements of density and energy are decoupled. 

We do have for the modes associated to the density, n{r,t), the characteristic equation 

iu(iu + ej^) - iQ ■ (siiQ + = 0 , (97) 

" \ m*/ 

or 

— iujOj^ — BiQ^ = 0 , (98) 

after neglecting the term with iQ ■ E (we recall that the electric held is constant and then 
its divergence is null). Solution of Eq. (98) provides us with the two roots, 

± i , (99) 

which can be rewritten as 

= 'fl! ± - 1 , ( 100 ) 

From this Eq. (100) we can characterize two types of movement: 

1 . 4BiQ‘^6]^ < 1, an overdamped regime, when uj± are purely imaginary, and in the limit 
AB^Q^el « 1, we can write 

<^±(Q) = ± (1 - . ( 101 ) 

and then 

^^+(Q) ^(102) 

= (103) 

where = Bidj^ is a diffusion coefficient. Hence, the hydrodynamic movement for Q 
sufficiently small is of the diffusive type. 

2. 4BiQ‘^6‘j^ > 1; then u±{Q) have an oscillating part and a relaxation time For 
ABiQ'^O'j^ 3> 1, we obtain 

UJ±{Q) - ± 1 (104) 

with having dimension of velocity. Hence, the hydrodynamic movement for 

sufficiently large Q is of the type of a damped wave, where v is the group velocity of the 
wave, and the dispersion spectrum is linear in the wavenumber. 
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It can be noticed that for any fluid a transition from one regime to the other (diffusion 
and damped wave) follows at a cut-off Qco given by . Movements well 

characterized by small wavenumbers {Q < Qco), are well described in MHT of order zero 
(the classical-Onsagerian one), i.e. by a Fick diffusion equation. Movements characterized 
by wavenumbers Q > Qco, are well describe in a MHT of order 1, implying in a damped 
wave equation (Maxwell-Cattaneo equation). This is up to a second cut off wavenumber, 
say Qi 2 , requiring for movements involving Q > Q 12 to go over a description in MHT of 
order 2 [38]. 

On the other hand, for the modes associated to the thermal motion we do have 


{iu -|- + ^if^) T C2Q^ — 0 , 

(105) 

which can be written as 

= 0 , 

(106) 

where 

= + ■ 

(107) 

The roots of Eq. (106) are 

^±{Q) = ± ^\Jc2Q^ + , 

(108) 

or 

^±{Q) - ^ 1 , 

(109) 

where 

A = 4C2T^Q^ + . 

(110) 


Quite similarly to the case of charge motion we just considered, we can evidence two 
regimes, namely 

• for M < 1, a purely diffusive regime, 

• for M > 1, a damped wave regime, 

with a cut-off wavenumber Qco dehning the frontier between both given by 

QL = (l-4f,-V,^)/402r^ (111) 

for values of Q < Qco there follows diffusive motion, and for Q > Qco damped wave motion. 
This is valid for any fluid, and, for example, can be visually observed in experiments on 
thermal stereolithography (or infrared laser induced rapid phototyping) [45]. 
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VI. CHARGE MOTION: ELECTRIC CONDUCTIVITY 

In Eq. (93), in the steady state and taking as null 634 and meaning that we disregard 
thermo-striction and diffusion effects, after integration in space we do obtain that 

el„. = aE, ( 112 ) 

where 

a = , (113) 

is the usual Sommerfeld-Drude expression for the electric conductivity; el„g is the electric 
current density and Eq. (112) is Ohm law. 

On the other hand, looking for the space dependence of the current, after differentiating 
on time Eq. (93) and using Eq. (92) there follows that 

n2 1 n 

eE 

—V-I„, (r, t) + Ge(r, t) , (114) 

m* 

resembling a Maxwell-Cattaneo-like equation and the telegraphist equation of electrody¬ 
namics, and Ge(r,t) is the contribution arising out of the term $ in the evolution equation 
for the density, namely 

G.(r,«)= (Bp''V + ^5)4(r,«). (115) 

V 

In the steady state and assuming isotropy such that after multiplying by 

Eq. (114) becomes 

ReV^lnXr) + SeEV ■ I„,(r) + 0/„,Ge(r) = 0 , (116) 

with Re = Oi^^Bie and Se = eOi^Jml. 

The general solution of Eq. (114) is a sum of solutions of the associated homogeneous 
one (obtained for Ge(r) = 0) and a particular solution with Ge(r) ^ 0. 

Consider Erst the solution of the homogeneous equation. Taking 2 : as the direction along 
the axis of the cylinder, and the electric held E parallel to it, and in cylindrical coordinates 
neglecting the dependence on the angle 9, and introducing a separation in variables r and 
2 , we have, after writing 

h{r,z) = /W(r) X lf\z ), 
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that 


1 d 


iPir) r 9r 




+ 


lP(z) 


IP {z)+^i-ip{z) 


dz'^ 


Rp dz 


= 0 , 


whose solution [46] is a sum in 7 of terms like 

h{r, z) = Ml, r) [A,e‘+'">* + , 


(117) 


(118) 


where 7 is a real number to be determined by the use of boundary conditions, Jo is Bessel 
function of the second kind, and ^±( 7 ) are the roots of 

k‘^ + ^^k + p = 0, (119) 

itg 


i.e.. 


k±p) 



0 . 


( 120 ) 


Notice that for 7 ^ ::§> SpS/Rp, k± i'y and the 7 -dependence of the solutions is oscil¬ 
latory. For 7 ^ SpS/Rp the two roots become k+ = 0 and k- = —SpE/Re^ hence the 
7 -dependence is “overdamped” for k_ and independent of 7 for kj^. For 2 : 3> Rp/SpS > 7 
and 7 r 1 the solution become independent of and r. Then, making = aS one 
recovers Eq. (112). 

In general a non-uniform current distribution may follow from the space-dependent ef¬ 
fects present in the term G of Eqs. (114) and (115). These are, as already noticed, the 
distribution of impurities, presence of imperfections, influence of weldings, and in the case 
of nanometric dimensions the question of boundary conditions in the presence of rugous 
walls with fractal on average topography, leading to a complicate reflection of the carriers 
[47]. Then the coefficients and in Eq. (118) may have to be adjusted for the complete 
solution to satisfy given boundary conditions. These situations are quite difficult to deal 
with theoretically, comprising a case of the so called ’’hidden constraints” in systems with 
complex structure [48]. An analysis of this question, i.e., the presence of G and complex 
boundary conditions, together with a study of the transient regime shall be reported in a 
future communication. 

Finally, according to the results here presented, a priori it appears that the conductivity 
is weakly dependent on the radius of the cylinder, but it is limited: The hydrodynamic 
treatment we have presented involve the motion of the average of a number of particles in a 
volume element, say (Rr, around position r. Hence, the results can not be extrapolated to 
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systems with very short nanometer dimensions, that is, involving lengths comprising a few 
lattice parameters: considering a, say, 5 A lattice parameter it can be suggested that the 
results are valid only for lengths larger than 50 to 100 nm. For smaller distances the motion 
would be greatly constrained and, as a rule, with the conductivity becoming much smaller 
than the one in bulk. 

VII. HEAT MOTION OF CARRIERS AND PHONONS 

The subject has been dealt with and reported in Refs. [31], and here we summarize the 
results for the sake of completeness of the topic. We consider Eqs. (94) and (95) for the 
carriers’ density of energy and its hrst flux, i.e., we take in direct space 

^he(r, t) + V ■ RAr, t) = nocE ■ t) - 0^^he{r, t) , ( 121 ) 

kl,.(r, t) = -Cf'Vfc.(r, i) - I»,(r, i) + t) , (122) 

where we have neglected thermo-electric effects, that is, we have taken &A, 6 I 3 , and as 
null. In the steady state they become 

V ■ RAr) = UoeE ■ RA^) - , (123) 

C'IVhe(r) = BIEKir) - , (124) 

where we have taken the tensors and C*! as scalars. 

From Eq. (122) we do have for the the heat current that 

RAr) = -C|0,,VhAr) + B^,ej„EK{v ), (125) 

and taking into account that the density of energy he can be written as 

heir) ~ ^nokBT*ir) + , (126) 

being composed of the thermal energy characterized by the nonequilibrium quasitemperature 
T*, and of the kinetic energy involving the drift-carrier velocity Ve- The latter as a general 
rule is smaller than the thermal energy [49] and disregarding it we can write 

I/ie(r) = -XeVTAr) + L,,E, (127) 
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can be interpreted as the carriers thermal conductivity which is space independent, and 



(129) 


can be considered as the carriers’ thermo-electric coefficient. Moreover 




(130) 


and 



(131) 


On the other hand, the thermal transport by phonons has been considered elsewhere. 


however in intrinsic semiconductors [31]. In doped semiconductors the influence of the 
electric held on the distribution of phonons is presented in Ref. [43], where it is shown the 


presence of a kind of resonance (overheating of certain reduced number of phonon modes in 
an off-center region of the Brillouin zone), which is not particularly relevant, arising out of 
the process of drifting electron excitation [44]. Therefore, we can state that the phonons’ 
thermal conductivity is very weakly affected by the presence of the electric held. 

The phonons’ thermal conductivity is strongly affected by the value of the radius of the 
cylinder in the nanometer domain. In Fig. 1 it is shown such dependence. Parameter £, 


o 

y 



0 5 10 15 20 40 60 80 100 


R/i 


FIG. 1: Dependence on the scaled wire radius, R/i, of the scaled thermal conductivity, k/kq, for 
several values of the reflection coefficient A and = s^OhOj^, after Ref. [31]. 








with dimensions of length, is a characteristic length, with = s'^OhOi ^, that is, in a Debye 
model, its square is given by the square of the sound velocity times the product of Maxwell 
times associated to the phonons energy density and energy flux. Figure 1 tells us that 
there follows a drastic reduction in thermal conductivity for R/t below the value 10, and 
becoming orders of magnitude smaller for R/^ < 1. We may then state that in the range of 
values oi RjI there exists a threshold below which the sample size (radius of the cylinder 
in units of i') leads to a notable reduction of the thermal conductivity, and large increase 
of the figure-of-merit in thermo-electric engineering. Figure 1 provides information on the 
influence of the reflection effect at the side boundaries: as expected with increasing reflection 
coefficient A there follows an increase in thermal conductivity. It must be noticed that we 
have considered normal reflection at an smooth surface, but the surface is always rugous 
with characteristics fractal on average [39] what affects the reflection processes. 

Taking into account the experimental data reported by D. Li et ah (Fig. 1(a) in Ref. [50], 
where it is shown the measured thermal conductivity of silicon in terms of the temperature) 
in samples of Si nanowires with different diameters (diameters of 22; 37; 56; and 115 nm), 
we consider those at 300 Kelvin, what is shown in Fig. 2. If we admit that for all the four 
samples kq is approximately the same and of the order of the thermal conductivity in bulk, 
namely kq — 148 (W/K.m) [51], we can obtain the values of k/kq given in Table I (third 
column), and from Fig. 1 (for M = 0, i.e., no reflection at the lateral borders: Couette-like 



FIG. 2: Measured thermal conductivity k of wires of Si in terms of the radius of the wire R, at 
300 Kelvin; experimental results (□) from Ref. [50]; after Ref. [31]. 
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flow) we can evaluate that, roughly, the corresponding values oi RjI are those given in the 
fourth column, and from them we can estimate the values of i shown in the hfth column. 
Considering as similar the Maxwell times for energy and its flux, which are equal in a Debye 
model, that is. Oh = Oj = 9, we get that y/Si/s ~ 6, and taking an average sound velocity 
of 8433 m/s, we obtain the values for the Maxwell time displayed in column 6 of Table I. 
The experimental data (open square dots) in Fig. 2 are contained in the curve (full line) 
adjusted by the second order polynomial 

K ~ -O.OUR^ + 1.65R - 9.31, 

for i? > 10 nm. The traced line for i? < 10 nm is an intuitive extrapolated indication, given 
by: K~ 0.045i?2 + o.l22i?. 


TABLE I: Results for Si 


R (nm) 

K (W/K.m) 

kIkq 

R/e 

i (nm) 

e (ps) 

11.0 

6.76 

0.046 

0.626 

17.57 

3.61 

18.5 

17.57 

0.119 

1.063 

17.40 

3.57 

28.0 

25.68 

0.173 

1.339 

20.91 

4.29 

57.5 

40.54 

0.274 

1.850 

31.08 

6.38 



FIG. 3: Measured thermal conductivity k of wires of GaN in terms of the radius of the wire R, at 
300 Kelvin; experimental results (Q) from Ref. [52]; after Ref. [32]. 
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From the experimental data reported by C. Guthy et al. (Fig. 2(a) in Ref. [52], where it 
is shown the measured thermal conductivity of GaN in terms of the temperature) in GaN 
nanowires with different diameters (diameters of 100; 126; 160; and 181 nm), we consider 
those at 300 Kelvin, what is shown in Fig. 3. Using the value of the thermal conductivity in 
bulk for GaN, namely Kq ~ 210 (W/K.m) [53,54] and taking an average (in this hexagonal 
crystal) sound velocity of 5170 m/s [52], we obtain, similarly to Table I, the values shown in 
Table II. The experimental values (open circular dots) in Fig. 3 are contained in the curve 
(full line) adjusted by the second order polynomial 

K ~ -0.0002R2 + 0.18i? + 5.7, 


for i? > 50 nm. 


TABLE II: Results for GaN 


R (nm) 

K (W/K.m) 

k/ko 

R/l 

^ (nm) 

9 (ps) 

50.0 

14.1 

0.067 

0.769 

65.0 

21.7 

63.0 

15.6 

0.074 

0.813 

77.5 

26.0 

80.0 

18.9 

0.090 

0.905 

83.4 

27.9 

90.5 

19.7 

0.093 

0.922 

98.1 

32.9 


It must also be noticed the important point that the characteristic length £ and Maxwell 
times depend on R and on the nonequilibrium thermodynamic state of the system. This is 
so because of their dependence on R, which determines the frequencies Unq^ and of the sum 
over nqz- Figure 4 shows the dependence on the wire radius R of the characteristic length 
^ and Maxwell time 9 for GaN and Si nanowires. The linear expressions that relate the 
characteristic length and Maxwell time with the radius R are indicated within the hgures. 

Goncerning the so-called hgure of merit, Z, which is a number that allows for obtaining 
a useful insight for optimizing design parameters, is constructed by choosing the parameters 
that are most centrally vital to a design solution. For the case of thermo-electric devices is 
used [3] 

Z = SealK, (132) 

where Se is Seebeck coefficient. If we consider Seebeck effect and the electric conductivity 
as nearly independent on size, and the phonon thermal conductivity as the relevant one, we 
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FIG. 4: Dependence on the wire radius R of the (a) characteristic length i and (b) Maxwell time 
9, for GaN and Si nanowires, after Ref. [31]. 


can see that, according to the resnlts in hgnres 1 to 3, the hgnre of merit Z of Eq. (132) 
greatly increases in qnantnm wires with radins in the interval of 10 to 90 nm. 


VIII. CONCLUDING REMARKS 

We have presented an extended theory of the Mesoscopic Hydro-Thermodynamics of 
phonons and carriers in n-doped direct gap polar semiconductors in the presence of electric 
helds. MHT, also referred to as Higher-Order Generalized Hydrodynamics, extends standard 
(or Onsagerian) hydrodynamics allowing to incorporate hydrodynamic motion not restricted 
to smooth in space and time characteristics (i.e., including intermediate to short wavelengths 
and intermediate to high frequencies). It consists in deriving a set of coupled hydrodynamic 
equations for the densities of quasi-particles (carriers and phonons) and of energy and their 
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fluxes of all orders. This has been done in Section III. 

The matter has been illustrated resorting to a MHT of order 1 for carriers and phonons, 
which is a contracted description in terms of their densities, energies, and the vectorial 
fluxes (electric current and heat cnrrent) of both. Criteria for performing snch contraction 
are discussed in Ref. [38]. 

The corresponding four hydrodynamic eqnations are coupled together, but if we disregard 
the cross-contributions associated to thermo-electric effects, there follows the separate sets 
of two equations for the motion of charges and two equations for the motion of energy. 
These are the basic Eqs. (69) and (70), and Eqs. (71) and (72), respectively. It may be 
noticed that in these equations are present the qnite important generalizations of Maxwell 
time. We recall that the origin of Maxwell time goes back to the fnndamental article by 
J.C. Maxwell in 1867, on the dynamical theory of gases and liqnids [27], in the strain rate 
model there presented it is considered as representing the time during which the stresses are 
damped [28]. Section IV is closed with an interpretation of the several contributions to the 
hydrodynamic equations. 

In Section V the hydrodynamic modes in this MHT of order 1 are derived. They allow 
to characterize the two regimes that are covered by it, namely, a diffusive motion at low 
wavenumbers and a damped wave motion at intermediate wavennmbers. In the first case the 
motion is governed by a typical diffusion eqnation (Fick’s and Fourier’s type respectively), 
and in the second by a Maxwell-Cattaneo-like equation. A cut-off wavenumber, Qco of Eq. 
(100), dehnes the frontier between the two types of regimes. 

Charge motion and characterization of the electric conductivity in the steady state are an¬ 
alyzed in Section VI. It may be notice that at no too small nanometric sizes the conductivity 
is nearly constant and taking a Drnde-type expression. Minor space-dependent effects may 
result from the presence of the space-dependent distribution of impurities, imperfections, 
inflnence of weldings, and boundary conditions (which have a ruggedness of a fractal-on- 
average type). As noticed at the closing of Section VI, the results can not be extrapolated 
to wires with very short nanometer radius, say, below a few tenths of nanometers. 

Heat motion and characterization of the thermal conductivity in the steady state are 
analyzed in Section VHI. In the case of the carriers, as it happens with the electric conduc¬ 
tivity, the thermal conductivity is constant with an expression of the type of standard kinetic 
theory, and taking into acconnt the expression for the electric conductivity there follows a 



type of Wiederman-Franz law. In the case of the phonons, quite differently, there follows 
a strong space dependence affected by the value of the radius of the wire. There follows a 
drastic reduction in the thermal conductivity as the radius decreases, evidenced within this 
MHT of order 1, which is being suppressed if one resorts to standard hydrodynamics. This 
may be interpreted that as the radius decreases to the nanometric scale larger wavenumbers 
need be included for the proper description of the movement. A MHT of higher order than 
1, would be required for wires with radius of a few nanometers. 

Finally, as a consequence of those results we can draw the attention to the fact that the 
so-called hgure of merit in the engineering of thermo-electric devices would greatly increase 
following the decrease of the wire’s radius. 
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Appendix A: The Nonequilibrium Statistical Operator 


According to NESEF ([9,11-13,25] with a short overview given in Ref. [32]), the nonequi- 
librium statistical operator in terms of the basic nonequilibrium variables in sets (8) and (9) 
is given by 

T^e{t) = Qeit) X qb , (Al) 

where 

t 

ge{t) = exp I In g{t, - J > (A2) 

— OO 

with g{t, 0) being the auxiliary statistical operator (also called “instantaneous quasi- 
equilibrium operator”) and 

p(f', t' -t) = exp I - ^ Ek(t')^k(t' -t) - ^ i^kQ(t')hkQ(t' - t) 

k k,Q5>^0 

- , (A3) 

q.7 q.Q^o 

7 

where t' stands for the dependence on the time of the nonequilibrium thermodynamic vari¬ 
ables F’s and the dynamical microvariables, in Heisenberg representation, depend on {f — t). 
Moreover gs is the canonical distribution of the bath of acoustic phonons in equilibrium at 
temperature Tq, and 0(t) ensuring the normalization plays the role of the logarithm of a 
nonequilibrium partition function. 

We recall that the second term in the exponent in Eq. (A.2) accounts for historicity 
and irreversibility in the nonequilibrium state of the system. The quantity e is a positive 
inhnitesimal that goes to zero after the trace operation in the calculation of averages has 
been performed. We also recall that 

ge{t) = g{t,0) + g'^{t) , (A4) 

i.e., it has an additive composition property, with a contribution of the instantaneous quasi¬ 
equilibrium statistical operator plus the one of which contains the historicity and produces 
irreversible evolution. 
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Appendix B: NESEF-Kinetic Theory 


The NESEF-based Kinetic Theory of relaxation processes basically consists into taking 
the average over the nonequilibrium ensemble of Heisenberg (or Hamilton at the classical 
level) equations of motion of the dynamical operators for the observables, say, Aj{r), with 
j = 1,2,..., (a function over phase space in classical mechanics and Hermitian operator in 
quantum mechanics) under consideration, i.e. 

= ^Tr{Aj{r)ge{t) x qb} = Tr !^^[Aj{r), H]ge{t) x , (Bl) 

which is a manifestation of Ehrenfest Theorem. The practical handling of this NESEF- 
Kinetic Theory is described in Refs. [9,11-13] and mainly in [26]. The NESEF is a powerful 
formalism that provides an elegant, practical, and physically clear picture for describing 
irreversible processes, adequate to deal with a large class of experimental situations, as for 
example, in semiconductors far-from equilibrium, obtaining good agreement in comparisons 
with other theoretical and experimental results [55]. 

Here we briefly notice that the Markovian limit of the kinetic theory is of particular 
relevance as a result that, for a large class of problems, the interactions involved are weak 
and the use of this lowest order, second order in the interaction strengths, in the equations 
of motion constitutes an excellent approximation of good practical value. By means of a 
different approach, E. B. Davies [56] has shown that in fact the Markovian approach can be 
validated in the weak coupling (in the interaction) limit. 

Explicitly written, the Markovian equations in the kinetic theory are 

^Aj{r,t) = (B2) 

where, after it is introduced in the Hamiltonian the separation H = Hq + H', where Hq 
stands for the kinetic energy and H' contains the interaction potential energies present in 
Eq. (Bl), we have that 

4 °^ ^ ^ ^ 
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and t) = t) + with 


_ t(2) 


( 2 ), 


= 77 


ihy 


^i'edt'-dTr I [H'(t' - t)o, [H', Mr)]] g{t, 0) x , (B5) 


= ;^ X] [ *^Tr |[if',ifc(r)]^(t,0) X 


SAk{r,t) 


(B6) 


where g is the auxiliary statistical operator of Eq. (A.3) and g^ the equilibrium statistical 
distribution of the thermal bath, and we recall that and M\ which in Mori’s terminology 
[57] are called the precession and force terms, are related to the non-dissipative part of the 
motion, while dissipative effects are accounted for in Jj which can be called scattering 
integrals. Subindex nought indicates evolution in the interaction representation, <5 indicates 
functional differentiation [44], 


Appendix C: Summary of Heims-Jaynes Procedure 


Given an statistical operator of the form 




A+B 


where 


Z = Tr{e^+^} , 

ensures its normalization, and introducing 


^?o — 


Tr{e^} ’ 

according to Heims-Jaynes, given an any operator 0 it follows that 

OO 

Tr{eq = (§)„ + 5](<5„(e - (§>„)), 

where 


n=l 


(0)o = Tr{0f?o}, 


with 


n—1 


Qn Sn ^ ^j{Qn}oSn—k j 


k=l 


(Cl) 

(C2) 

(C3) 

(C4) 

(C5) 

(C6) 
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for n > 2 , and Qo = 1 and Qi = Si, 

TDn 

Sn = ^, So = l. (C7) 

n\ 

Equation (C4) consists of the average value of 0 with (that is, only depending on 
A) plus a contribution in the form of a series expansion in powers of B. In a hrst-order 
approximation we do have that 

Tr{0^?} ^ (0)o + Tr{5(0 - (0)o)^?o} • (C8) 

In Section III we have used that 

0) = . (C9) 

where 

A{t) = Fn^{t)ne + ■ In^ + 

Fhe{t)he + Fh^{t) ■ Ih^ + 

Fnp{t)np + F,T,p(t) ■ + 

Fhp{t)hp + Fh^{t) -Ihp , (CIO) 

that is, the homogenous part, Q 7 ^ 0, in the exponent of Eq. (48), and B(t) is the inho¬ 
mogeneous part, meaning the contributions with Q 7 ^ 0, and we have used the hrst-order 
(linear in B) approximation. 

In particular we had that 

^J(Q. i) = B'S(t)n,(Cl,t) + B2(i)A.(Q, i), (Cll) 

4hQ.«) = c2(i)n7Q.«) + c2(i)fc.(Q.«). (C12) 

with tensor B and C given in Eqs. (73) to (76). 

On the other hand, for the case of the phonons we do obtain that 


^'q(^) - (^)] [Fn (t) ' VqWq "2 Ffe(t) ■ /hUqVqWq] , (C13) 


where 


'^q(^) = 


1 




(C14) 


33 



that is, a first order Taylor expansion in F„ and F/i (linear approximation). 

Next, resorting to the use of the nonequilibrium equations of state that relate the four 
nonequilibrium thermodynamic variables to the four basic variables, it follows in hrst-order 


Heims-Jaynes expansion that 

n(Q,t) = Auit)(pniQ,t) + Ai 2 it)(phiQ,t ), (CIS) 

I„(Q,i) =43(i) ■F«(Q.i) +44(i) ■F).(Q,«). (C16) 

h(Q, t) = A2i(t)v9n(Q, t) + A22it)(phiQ, t) , (C17) 

I,(Q,t) = 41 (^) ■ F„(Q,t) +41(^) ■ F„(Q,t), (C18) 


where An, A^, 41) 41) ^ 21 , A 22 , A^l and 41 those of Eqs. (C24) to (C29) below, 
except for the replacement of Z 2 q(t) of Eq. (C13) by Fq(t) of Eq. (C14). 

In Eqs. (CIS) and (C17) the contributions in F„ and F/^ present in Eq. (C13) are null, 
whereas in Eqs. (C.16) and (C.18) are null the contributions in (pn and ph- Eqs. (CIS) 
to (C18) constitute a set of linear algebraic equations that can be inverted to obtain the 
four nonequilibrium thermodynamic variables pn, Ph, F„ and F/j, in terms of the basic 
hydrodynamic quantities, n, h, I„ and I^. 

The second-order fluxes are given by 

/f''(Q,«)=3g(iK(Q,«)+lMW*)»(Q.i). (C19) 

4'(Q.<)=3S(iK(Q,«)+44Ws)/.(Q.i). (C20) 

_[ 2 ] _[ 2 ] _[ 2 ] 

where Agg, A^^ and A^^ are those of Eqs. (C27), (C28) and (C29), except for the replacement 
of Z2q(t) of Eq. (C13) by Fq(t) of Eq. (C14). 

On the other hand, introducing the concept of nonequilibrium temperature, better called 
quasitemperature T*{r,t) in the form 

= (C21) 

we can obtain an evolution equation for it starting with the evolution equation for the 
energy in the form of the hyperbolic Maxwell-Cattaneo equation, from which together with 
the nonequilibrium thermodynamic equation of state, Eq. (C17), we have that 

^(ha;q)2z7q(t)[l + Z7q(t)] 


d‘^Ph{r,t) 

dt^ 


+ 4 


, dph{r,t) 


+ 


h{r,t) 
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= V 


3g(o4H§-:42(o 


Viph{r,t) + 


41(«) 


^11 (^) 

and, after introducing the heat capacity 


Vn(r,() + »r'ir“(r,(), 


(C22) 


Cv{t) = kBY^ 

where Tq is the temperature in equilibrium in this linear treatment, and the quantities Aij 
and Aij are. 


^n{t) = ^Z/q(t)[l + Z/q(t)], 

q 

(C24) 

Auit) = A 2 l{t) = ^z/q(t)[l + z/q(t)]ha;q, 

q 

(C25) 

dl22(t) = ^Z/q(^)[l + Z/q(^)](h4Jq)^ 

q 

(C26) 

^y(^) ~ ^q(^)[l + ^q(^)] [VqWqVqCUq] , 

q 

(C27) 

= 41 W = '^qWll + ^^qWllVqWqVqWqj/h^q, 

q 

(C28) 

41 (^) = Yl + ^q(^)][VqWqVqa;q](h4Jq)^ , 

q 

(C29) 

^12(^) = ^ll(^)^22(^) ~ ^12(^)^12(^) 5 

(C30) 

A34(t) = 41 w © 41(^) - 41(^) © 41 w • 

(C31) 


In these expressions, [VqCaqVqWq] denotes the second order tensor with components 

duj/dqiduj/dqj, while = '^FijGji^ with © standing for full contracted description. 

ij 


35 



[1] C.G. Rodrigues, A.A.P. Silva, C.A.B. Silva, A.R. Vasconcellos, J.G. Ramos, R. Luzzi, Braz. 
J. Phys. 40(1), 63 (2010). 

[2] D. Jou, J. Gasas-Vazquez, M. Criado-Sancho, Thermodynamics of Fluids Under Flow 
(Springer, Berlin, Germany, 2001). 

[3] D.M. Rowe, Thermoelectric Handbook: Macro to Nano (Taylor and Francis, Boca Raton, 
USA, 2006). 

[4] Yu L. Klimontovich, Statistical Theory of Open Systems: A Unified Approach to Kinetic De¬ 
scription of Processes in Active Systems, vol. 1, Kluwer Academic, Dordrecht, The Netherlands 
(1995). 

[5] S.K. Belyaev, Phys. Uspekhi 38, 287 (1995). 

[6] D.N. Zubarev, V.G. Morozov, I.P. Omelyan, M.V. Tokarchuk, Theor. Math. Phys. 96, 997 
(1994). 

[7] B.J. Adler, D.J. Tildesley, Computer Simulation in Liquids, Oxford Univ. Press, Oxford, UK 
(1987). 

[8] M.H. Kalos and P.A. Whitlock, Monte Carlo Methods, (Wiley Interscience, New York, USA, 
2007). 

[9] D.N. Zubarev, Nonequilibrium Statistical Thermodynamics, Plenum-Consultants Bureau, New 
York, USA, 1974. 

[10] J.P. Dougherty, Phil. Trans. Roy. Soc. (London) A 346, 259 (1994). 

[11] R. Luzzi, A.R. Vasconcellos, J.G. Ramos, Predictive Statistical Mechanics: a Non-Equilibrium 
Ensemble Formalism, Kluwer Academic, Dordrecht, The Netherlands (2002); and Springer e- 
Books Archive. 

[12] R. Luzzi, A.R. Vasconcellos, J.G. Ramos, Rivista Nuovo Cimento 29(2), 1-85 (2006). 

[13] D.N. Zubarev, V. Morozov, G. Ropke, Statistical Mechanics of Non Equilibrium Processes, 
vols. 1 and 2, Academie Verlag-Wiley VCH, Berlin, Germany (1996). 

[14] L. Sklar, Physics and Chance: Philosophical Issues in The Foundations of Statistical Mechan¬ 
ics, Cambridge Univ. Press, Cambridge, UK (1993). 

[15] C.G. Rodrigues, A.R. Vasconcellos, R. Luzzi, Phys. Stat. Sol. (b) 246(2), 417 (2009); ibid. 
Phys. Stat. Sol. (b) 216, 35 (1999); ibid. Braz. J. Phys. 32(2A), 439 (2002); ibid. Braz. J. 


36 



Phys. 36(2A), 255 (2006); ibid. Appl. Phys. Lett. 82, 2455 (2003). 

[16] C.G. Rodrigues, A.R. Vasconcellos, R. Luzzi, J. Appl. Phys. 102, 073714 (2007); ibid. J. 
Appl. Phys. 90, 1879 (2001); ibid. Solid State Commun. 140, 135 (2006); ibid. Eur. Phys. J. 
B 72(1), 67 (2009); ibid. J. Phys. D - Appl. Phys. 38, 3584 (2005). 

[17] Chapter 6 in the book of Ref. [11]. 

[18] H.J. Kreuzer, Nonequilibrium Thermodynamics and its Statistical Foundations, Clarendon, 
Oxford, UK (1981). 

[19] H.J.B. Casimir, Rev. Mod. Phys. 17, 343 (1945). 

[20] D. Jou, J. Casas Vazquez, C. Lebon, Extended Irreversible Thermodynamics, Fourth Edition, 
Springer, Berlin, Cermany (2010). 

[21] I. Muller, T. Ruggeri, Extended Thermodynamics, Springer, Berlin, Cermany (1993). 

[22] D. Jou, J. Casas-Vazquez, C. Lebon, Rep. Prog. Phys. 51, 1105 (1998); ibid. 62, 1035 (1999). 

[23] J.P. Boon, S. Yip, Molecular Hydrodynamics, McCraw-Hill, New York, USA, (1980); reprinted 
by Dover, New York, USA (1991). 

[24] T. Dedeurwaerdere, J. Casas-Vazquez, D. Jou, C. Lebon, Phys. Rev. E 53(1), 498 (1996). 

[25] A.I. Akhiezer, S. V. Peletminskii, Methods of Statistical Physics, Pergamon, Oxford, UK 
(1981). 

[26] L. Lauck, A.R. Vasconcellos, R. Luzzi, Physica A 168, 789 (1990). 

[27] J.C. Maxwell, Phil. Trans. Roy. Soc. (London) 157, 49 (1867). 

[28] L.D. Landau, E.M. Lifshitz, Theory of Elasticity (Pergamon, London, UK, 1959). 

[29] C.A.B. Silva, J.C. Ramos, A.R. Vasconcellos, R. Luzzi, J. Stat. Phys. 112, 2692 (2011). 

[30] C.A.B. Silva, J.C. Ramos, A.R. Vasconcellos, R. Luzzi, Mesoscopic Hydrothermodynamics: 
Foundations Within a Nonequilibrium Statistical Ensemble Eormalism, arXiv:1210.7280. 

[31] C.G. Rodrigues, A.R. Vasconcellos, R. Luzzi, Eur. Phys. J. B 86, 200 (2013); C.G. Rodrigues, 
A.R. Vasconcellos, R. Luzzi, Physica E 60, 50 (2014). 

[32] F.S. Vannucchi, A.R. Vasconcellos, R. Luzzi, Int. J. Modern Phys. B 23, 5283 (2009). 

[33] U. Fano, Rev. Mod. Phys. 29, 74 (1957). 

[34] N.N. Bogoliubov, Lectures in Quantum Statistics I, Gordon and Breach, New York, USA 
(1967). 

[35] J.R. Klauder, B.S. Skagerstam, Coherent States, World Scientific, Singapore (1984). 

[36] N. Hugenholtz, Application of Field-Theoretical Methods to Many-Boson Systems, in 1962 


37 


Cargese Lectures on Theoretieal Physies, M. Levy Ed., Benjamin, New York, USA (1963). 

[37] R. Balian, Y. Alhassed, H. Reinhardt, Rhys. Rep. 131, 1 (1986). 

[38] J.G. Ramos, A.R. Vasconcellos, R. Luzzi, J. Chem. Phys. 112, 2692 (2000). 

[39] F. Family, T. Vizsek, Eds., Dynamic of Fractal Surfaces (World Scientific, Singapore, 1994). 

[40] C.G. Rodrigues, V.N. Freire, A.R. Vasconcellos, R. Luzzi, Mat. Res. 6(1), 1 (2002). 

[41] G.G. Rodrigues, A.R. Vasconcellos, R. Luzzi, V.N. Freire, J. Phys.: Condens. Matter 19, 
346214 (2007). 

[42] S.P. Heims, E.T. Jaynes, Rev. Mod. Phys. 34(2), 143 (1962) (see Appendix B in page 164). 

[43] G.G. Rodrigues, A.R. Vasconcellos, R. Luzzi, J. Appl. Phys. 113, 113701 (2013) 

[44] G.G. Rodrigues, A.R. Vasconcellos, R. Luzzi, J. Appl. Phys. 108, 033716 (2010). 

[45] J.G. Ramos, A.R. Vasconcellos, R. Luzzi, Int. J. Quantum Ghem. 65, 277 (1997). 

[46] R. Courant, D. Hilbert, Methods of Mathematical Physics (Wiley-Interscience, New York, 
USA, 1953). 

[47] J.M. Ziman, Electrons and Phonons (Glaredon, Oxford, UK, 1960). 

[48] R. Luzzi, A.R. Vasconcellos, J.G. Ramos, Rivista Nuovo Cimento 30(3), 95 (2007). 

[49] G.G. Rodrigues, A.R. Vasconcellos, R. Luzzi, J. Appl. Phys. 99, 073701 (2006). 

[50] D. Li, Y. Wu, P. Kim, L. Shi, P. Yang, A. Majumdar, Appl. Phys. Lett. 83, 2934 (2003). 

[51] M. Asheghi, Y.K. Leung, S.S. Wong, K.E. Goodson, Appl. Phys. Lett. 71, 1798 (1997). 

[52] G. Guthy, C-Y Nam, J.E. Fischer, J. Appl. Phys. 103, 64319 (2008). 

[53] A. Jezowskia, P. Stachowiaka, T. Suskib, S. Krukowskib, M. Bockowskib, 1. Grzegoryb, B. 
Danilchenkoc, Physica B; Gondensed Matter 329-333 Part 2, 1531 (2003). 

[54] A. Jezowskia, B.A. Danilchenkoc, M. Bockowskib, 1. Grzegoryb, S. Krukowskib, T. Suskib, 
T. Paszkiewicz, Solid State Gommun. 128, 69 (2003). 

[55] G.G. Rodrigues, A.R. Vasconcellos, R. Luzzi, V.N. Freire, J. Appl. Phys. 98, 043703 (2005); 
G.G. Rodrigues, J.R.L. Fernandez, J.R. Leite, V.A. Chitta, V.N. Freire, A.R. Vasconcellos, 
R. Luzzi, 95, 4914 (2004). 

[56] E.B. Davies, Communications in Mathematical Physics 39(2), 91 (1994). 

[57] H. Mori, Progress of Theoretical Physics (Japan) 33, 423 (1965). 


38 



